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Abstract — A Goody random band model with a Voigt line profile is used to calculate the band 
absorption of CFCL at various pressures at room and stratospheric (216°K) temperatures. Absorption 
coefficients and line spacings are computed. 

1. INTRODUCTION 

The vibration-rotation bands of freon II (CFCK) have been studied in the past. 1 ' 5 This heavy 
symmetric top molecule with very small rotational constants has a multitude of rotational lines 
in the bands near 10 /im, which should explain the lack of any discernable rotational structure 
at modest resolution. Thus, the spectral absorption coefficients for trichloro-fluoromethane 
should be a slowly varying function of wavelength, making it an ideal candidate for a band 
model analysis. 


2. THEORETICAL CONSIDERATIONS 

The general theory of the random band model was first developed by N. Mayer 7 and later by 

Goody 6 and is only briefly described here. The absorption lines are assumed to be randomly 

spaced with their intensities distributed according to a Poisson distribution. The wings of the 

lines in the center of the spectral interval being modeled must have a negligible contribution 

outside that interval. The general form of the random band model is 


T = exp 


1 r kpl8f v 

8 J-x I + kpldfy 



( 1 ) 


where p is the density (amagat), / is the length (cm), k is the mean absorption coefficient (cm -1 
amagar 1 ), 8 is the mean line spacing (cm -1 ), and f v is a line shape function. Over the interval 
that defines T, fc, and 8 are constant. The integration over the Poisson distribution has already 
been performed. 

For the widest range of utility, the Voigt line shape is the most realistic profile to use for /„ . 
The appropriate absorption profile can be written (see Penner 7 ) as 


= # 2 ) 
h «dV(tt) 


V(x,y) = 


V(fn 2) y 

a D \/(7T)TT 


f x exp (-f 2 * * ) 

J_.y 2 + U-r ) 2 


d t 


( 2 ) 


where 

'-© Vfln2) - 

with a L and a D the Lorentz and Doppler halfwidths at halfmaximum, and v 0 is the frequency of 
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Ffierlussi et aT® The ^ ^ calculation ° f the V oigt Profile is available from 

Pierlussi et al. The integral in the general equation, Eq. (1), can be expressed as 


T = exp - 


Ikpl 

V(tt) 


/; 


V(x,y)dx 




(3) 


The Lorentz halfwidth can be expressed as 


a L ~ 0>„ a P l , + Otf°Pf, 

fcm-'/ a ^Hnr r t e hi he h Par K al Pre f K 6 ! (3tm) 3nd a “°’ are ,he pressure broadening coefficients 
(cm /atm) for the absorber and the foreign gas respectively. 

. ^ 3) was evalu ated by Romberg’s numerical integration technique as described by 

oor. The band model was fitted to the laboratory transmission data of Nanes et al 5 at 
constant temperature using a least squares fit to a non-linear function with a linearization of the 
fitting function as described in Bevington. 10 Using this method it is possible to determine the 
following parameters for a g.ven frequency interval, Av. k, mean absorption coefficient 

L- 1/ ™ 8 " I ’ TT l,ne SP u Cm8 <Cm ' ); " a °’ L ° rentZ self_broat * en > n g coefficient 
(cm / atm), a; , Lorentz foreign gas broadening coefficient (cm~ ! /atm) 

ca. bTexp°ess e d » “ ° bSerVa ' i0nS Can * meas " red * chi square rest which 


P-jJ-Stfl-fUiir 


(5) 


where T, ,s the measured transmission level. T(x,) is the calculated level, and N, is the number 
of degrees of freedom (number of spectra - number of parameters being determined). 

,, US ' n . 8 ^ q . (5) * * as calculated for given values of k, 8 , and a L °. These parameters were 
n adjusted until * was minimized. We have selected as our measure of uncertainty in the 

minimumy^ tvl^Th r $ "* UnCCrtain,y that co ^sponds to an increase of the 

,11 ,i!! Tk* by 100% ’ Tha IS ’ ‘ f WC Change one P arameter a m by an amount Aa m and optimize 
old value- ^ ParametCrS f ° r 3 new minimum then the new value of will be twice the 


X'(a m + Aa m ) = 2 x 2 {a m ). 

The utility of the band model can be extended by applying the following 

= f K dy — 2 (^)fAi >i. 

J i 


Thus, once the absorption coefficients have been determined, 
derived from those parameters of the model. 


the band intensity, S v °, 


may be 


3. RESULTS 

The average transmission for a given interval was calculated from the equivalent width 
by Nanes et aL that interval. This transmission value along with the experimental 
i l0ns o pressure, path length, and temperature formed the data sets for the modeling 
program. The computed results are presented in Table 1 and displayed in Fig. 1 for both the 

n c ^ and thC ,0W T 6±2 ° K) tCmperatUre Spectra ' The four ba " d complexest the 

\ a , r T° n Were m0deled at l0W reso,ution ( a broad interval covering an entire 
band) and at moderate resolution (intervals spanning 10 cm' 1 ). 

siJL™. ' ine Sh , aP r WCre Lorentzian ’ then Eq - (1) can be analytically integrated and a 
g cant amount of computation can be eliminated. Our least squares fitting algorithm failed 
to converge to any meaningful values when using the Lorentz line shape. Only when the Voigt 
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Fig. 1. The mean absorption coefficient, k, for freon 1! in the 770-1 120cm' 1 region at 294 and 2I6°K. 

profile was incorporated into the model was convergence achieved. The pressure domain of the 
data, 7xl0~ 6 to 3x|0 -4 atm, necessitated the use of the Voigt profile with the Doppler 
halfwidth of 6xl0" 4 cm -1 as the dominant term. Except for the 810-880cnr' region of the 
spectrum, where k is the largest, we were unable to determine a Lorentz halfwidth. For this 
region, however, we found 

0= [0.08 ±0.01 cm” '/atm at 294°K, 

L 10. 10±0.04cm '/atm at 216°K. 

These values were used for the Lorentz broadening coefficient in the band model in order to 
calculate the mean absorption coefficient and the mean line spacing for all other intervals. Over 
the pressure range of the available data, the model is insensitive to the Lorentz broadening 
coefficient. For the 8 10-880 cm' 1 region k changes by less than 6% as a L ° approaches zero. 
Because there was no significant difference detected when N 2 was added to the CFCIj, 5 we did 
not attempt to determine the nitrogen broadening coefficient. 

For a few of the intervals, the fitting routine was unable to converge upon a value for the 
mean line spacing. In these cases the value of 8 given in Table I is believed to represent an 
upper limit for the mean line spacing and therefore no uncertainty is given. For all intervals 
where absorption was observed, a value for the mean absorption coefficient, k, was determined. 
All the random errors in the laboratory data are reflected in the uncertainties (sigma) of the 
parameters listed in Table 1. If there are any systematic errors in the data, they must be 
factored into the results. 

The band intensities were determined using Eq. (7) based on the absorption coefficients for 
the 10cm“ intervals. The results are presented in Table 2. There is as much as a 6% difference 
between these intensities and those derived from the broad intervals. It was found, for the 
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Table 2. Absolute intensities, S° (cm , /cm-amagat),+ of the i.r. bands of CFCL (freon 1 1) in the 8-12 /im region. 


Interval 

Temperature 




(cm ^/cm-amagat) 


-1 








cm 

K 

Previous Work 

Band 

Model 



(ref.) 





294 

20.8 

+ 

3.6 

(5) 

24.0 

± 1.9 

760-810 

216 

11.2 

+ 

3.1 

(5) 

21.7 

± 7.2 



1688 

± 

49 

(1) 




294 

1669 

1876 

+ 

+ 

135 

281 

(2) 

(3) 

1734 

± 13 



1701 

+ 

48 

(5) 



810-880 

216 

1556 

± 

57 

(5) 

1641 

± 21 


294 

54.7 

+ 

2.7 

(5) 

55.5 

± 1.8 

910-960 

216 

49.2 

+ 

1.7 

(5) 

51.0 

± 1.2 



698 

♦ 

40 

(1) 




294 

781 

576 

i 

67 

(2) 

(4) 

650.0 

± 3.5 



645 

± 

15 

(5) 



1030-1120 

216 

624 

± 

17 

(5) 

638.1 

+ 3.2 


t 


etp 


is taken as equivalent to amagat 


810-880 cm" 1 region, that as the interval widths were narrowed, the intensity converged to the 
sum from the 10 cm 1 intervals. The intensity derived from the 10 cm" 1 intervals was well 
within the error limits for the summation based upon 20 cm" 1 intervals. Consequently, it was 
felt that the intensity derived from intervals narrower than 10 cm" 1 would be essentially the 
same as those reported in this work. This may be an artifact of replacing an integral with a finite 
sum in Eq. (7). 


4. DISCUSSION AND CONCLUSION 

The band model calculations accurately predict the laboratory observations, as evidenced by 
the low values of \ 2 in Table 1. Except for the very weak region, 760-810cm~\ at 216°K the 
agreement with the measured intensity values listed in Table 2 is excellent. The agreement can 
also be seen in Fig. 2 where the synthetic spectrum is superimposed upon an actual laboratory 
observation. It was determined that the placement of the intervals did not significantly effect 
the results for the integrated band intensities. 

The model predicts the measured transmission such that the residual |f obs - f cak |, is less 
than 0.005 in 90% of the cases and in no instance is the residual greater than 0.02. Having a 
model that is in excellent agreement with observations warrants a deeper examination of its 
paramefers. 

A comparison of the laboratory data and the band model calculations is displayed in Figs. 3 
and 4. It can be seen from these curves that as kpl increases, the value of the mean line spacing 
has more influence on transmission. For different regions of the 810-880 cm" 1 complex, when 
kpl = 1 the range of mean line spacings affects the transmission by as much as 10% (f is near 
50% at this point). We do not know of any measurements of the mean line spacing for CFCI3. 
However, Jennings 11 has measured freon 12 (CF 2 CI 2 ) with lO^cm" 1 resolution and has been 
unable to find any single isolated lines for this lighter molecule. His results are consistent with 
the very small values of the mean line spacing derived by this model. 

The other interesting parameter in the model is the Lorentz halfwidth, a L °. Despite a less 
than ideal set of laboratory conditions, we have derived a single value for both 294° and 216°K. 
It is generally accepted that a L ° is proportional to (T)~ m y with which our values are consistent, 
eg- 


a L °(294°K) 

a L 6 (216°K) 


0.8 ±0.3 


[V(2 16/294) = 0.86]. 
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Fig. 2. Spectrum of freon 1! in the 8-12 /xm region at 2 16 Q K with the synthetic spectrum derived from a 

band model superimposed. 



KpV 

Fig. 3. Comparison of laboratory data and band-model calculations for freon 11 in the 810-880 cm 1 region 

at 294° K. 


The similarity of transmission as a function of kpl for the 10 cm 1 intervals covering the 
band centers at 847 and 1085 cm 1 can be seen in Fig. 5. These two intervals exhibit the same 
behavior for both the room temperature and stratospheric temperature data. The small 
temperature dependence combined with the large absorption coefficients for these intervals 
should make them ideal candidates for stratospheric freon 1 1 abundance determination. 
Because these are 10 cm -1 intervals, low resolution experiments could be employed to take 
advantage of the results presented in this work. Band model parameters could be derived for 
intervals of higher resolution (0.25 cm - ’ maximum) with the current data, if required. 

This model compares directly with the laboratory measurements of Nanes et al . 5 over a 
range of 4000 kpl. With the extremely small mean line spacings found it seems that the band 
model could be applied to intervals of much higher resolution. It is hoped that his model will 
prove useful in atmospheric calculations. 
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Fig. 4. Comparison of laboratory data and band-model calculations for freon II in the 810-880 cm 1 region 

at 216°K. 



Fig. 5. Temperature dependence of band centers at 847 and 1085 cm 1 of freon 11. 
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